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Abstract 

The mixed model of the large N induced QCD, with Nj <C N flavors of heavy 
fermions in fundamental representation, is solved in the local limit. The Z/y sym- 
metry is broken spontaneously in the large ./V limit, evading the Elitzur "no-go" 
theorem. As a result of this symmetry breaking, there is the Bose condensate of 
the eigenvalues of the scalar field, proportional to This condensate leads to the 
mass unit, which goes to zero as fractional power of ^f, thus defining the new kind 
of the local limit of this lattice theory. There is a strong coupling region below this 
mass scale, which revives the hopes of induction of realistic QCD. 



1 Introduction 



In the previous paper [T!| we suggested the mixed model of induced QCD, with the local 
Zn symmetry broken by the extra spinor matter in the fundamental representation of the 
SU(N) gauge group. In the old model % P]| , built around the scalar field in adjoint 
representation of the gauge group, this symmetry led to the infinite string tension || ||, 
which was the main obstacle for the induction of QCD. 

Let us discuss some unusual features of induced theories, before going into technical 
details. 

It is not difficult to give an example of the model of induced QCD with broken sym- 
metry and correct perturbative behavior at intermediate scales.^ Just take large number 
Nf ^> N of flavors of the fermionic matter \&* , i — 1, . . . Nf in the fundamental represen- 
tation. Then, the effective QCD Lagrangean tr F 2 U would be induced with the factor 
Nf, which corresponds to small effective gauge coupling Ng 2 ~ jj- in the intermediate 
region. 

Unfortunately, such model would not be solvable. There would be about N 2 gauge 
invariant composite fields ^ X ^ 3 X at each cite x of the lattice after integrating out the gauge 
field. The effective Lagrangean for these fields would be only about TV 2 ( which is the 
number of the gauge degrees of freedom integrated out). So, the effective theory would be 
strongly coupled, even stronger than induced QCD. For the solvable model, this should 
be the other way around. The number of flavors should be much less, than the number 
of colors, so that the effective theory for the composite fields would be weakly coupled. 

On the other hand, we want the induced QCD coupling to scale as iV°, which requires 
about N 2 matter degrees of freedom. We could achieve this by adding some field in 
the adjoint representation. The unique property of the scalar field $ x with one flavor is 
that most of its N 2 degrees of freedom can be absorbed into the gauge field. Only TV 
eigenvalues of this matrix are left, the rest can be gauged away. One can go one step 
further, by introducing the eigenvalues density p x (4>) as collective field, by local change of 
variables ||. Then the effective Lagrangean is N 2 times greater than the number of the 
local degrees of freedom, so that the WKB expansion could be developed around some 
classical solution for p. So, the model could be solved for arbitrary dimension D of the 
lattice. 

Does this mean, that the model is trivial? Not necessarily. The classical solution for 
this effective field theory involves the quantum effects of original theory, since the majority 
of degrees of freedom is integrated out. Should we instead integrate out the matter field 
we would still get the induced gauge theory with correct counting of N. 

The question is, whether we are dealing with the correct phase of the gauge theory, 

1 There are two strong coupling regions in induced models, one at the lattice scales, and another at 
physical hadron scales, with the perturbative region in between. In the local limit, the whole range of 
scales expands, so that the lattice "confinement region" becomes unobservable. 
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in other words, are we able to reproduce the perturbative QCD at intermediate scales. 
This question is a dynamical one, there are no kinematical restrictions, which would 
immediately lead to the negative answer. It is crucial though, that there is the self- 
interaction of the scalar field. Without this interaction one gets only the lattice artifacts, 
as it was observed in numerical simulations as well, as in the analysis of the exact 
solution [U . 

With the self-interaction present, there is the following scenario of induced QCD |2|]. 
Let us assume we induced QCD and let us go backwards in spatial scales, from the large 
ones to the lattice spacing. The confinement and perturbative QCD regions are the same 
as usual, all the way to the scale when the scalar field revives. When we go from this region 
deeper into the ultraviolet direction, the effective scalar coupling grows, and destroys the 
asymptotic freedom.^. Some adjustment of parameters is needed, of course, otherwise we 
would not reach the exactly infinite gauge coupling at the lattice scales. Also, the bare 
scalar mass should be adjusted, otherwise the scalars would decouple at the lattice scales, 
rather than in the perturbative region. 

The density of eigenvalues of the scalar field reflects the evolution of the effective 
coupling as follows. The region of small (large) density correspond to strong (weak) 
effective coupling of the effective field theory, because the coupling is proportional to the 
ratio of the fluctuative part to the classical one (or, to put it another way, the ordered 
eigenvalues fluctuate only between neighbors, so, the denser they are spaced, the less they 
fluctuate). 

Since the scalar field has positive dimension of mass at D > 2, the scales of the 
eigenvalues are directly related to the mass scale. The induced QCD scenario would 
correspond to small density p((p) —>■ at small and large 0, with the large density region 
in between. The large strong coupling region is familiar. This is the string singularity, 
corresponding to the endpoint of the support of the eigenvalues. The double scaling limit 
of the matrix models of the string theory corresponds to the rescaling of the vicinity of 
this endpoint at large N. 

The small <fi strong coupling region is new. It corresponds to the infrared, rather than 
ultraviolet limit, where there are some massless excitations present. This corresponds to 
the confinement region in the induced QCD scenario. As for the large density in between, 
this must be the region of the perturbative QCD, if all goes according to the plan. 

The exact solution of the model, for arbitrary potential, to some extent agrees with 
this scenario. There are infinitely many fixed points, which correspond to different power 
laws p((j)) — > |0| a ,cos7ra = — 3 ^_ 2 for the density at the origin. The different values of a 

2 For this end, the screening of the gauge coupling by the interacting scalar field must overcome the 
usual antiscreening. In absence of scalar interaction, the screening is like 20 times smaller, than the 
antiscreening, as it is well known. Presumably, the strong attraction is needed, to put together 20 times 
more scalars. The attraction of scalars in the matrix models is a well-known phenomenon. Within the 
jf expansion it does not lead to instabilities, but beyond that it requires nonperturbative stabilization 
by higher terms. 
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correspond to the different relevant terms in the scalar potential. Namely, the term $ 2m+2 
contributes 2m+1 to the classical equation for density, which is relevant for a > 2m + 1. 
Thus, for 1 < a < 3 only the mass term is relevant, so that the mass should be adjusted. 
The quartic coupling shows up in the next level solution, with 3 < a < 5, etc. 

The simple dimensional counting, with $ ~ r 1 ^^ , shows, that the $ 2m + 2 terms are 
relevant in D < 2 + — dimensions. One, therefore, may speculate, that the fixed point 
with 2m + 1 < a < 2m + 3 is realized at 2 H — %-r < D < 2 + — dimensions. Thus, at 

m+l — m ' 

4 < D < oo the fixed point with 1 < a < 3 is realized, at 3 < D < 4 the one with 
3 < a < 5, etc. At D — > 2 + the power of the potential grows, so that it becomes 
nonpolynomial at D — 2. 

Formally, these solutions for the density apply for D = 1 as well, though, of course, 
the ultraviolet and the infrared regions are interchanged here, as the scalar field has 
the negative mass dimension. There, the gauge field can be gauged away, so that the 
model reduces to the C = 1 string theory with the well-known exact solution: p((p) = 



<2(E — U(4>)). In our solution, a = 2k + 1, which corresponds to the upside-down 
potential Uq — $ 4fe+2 ) with the chemical potential at the tip: E = Uq. f] 

The physical mass spectrum, corresponding to the plane wave fluctuations of the 

vacuum density, also satisfies the scaling law m 2 vhys oc e s , where e is the deviation of the 

1 

2 1 2(a-k) 

in |0. 



bare parameters from the critical values, and 5 = \ + \, 2 ™ k } ■ The details can be found 



These arguments applied to the scalar model. In the mixed model, with the perturba- 
tions, coming from the fundamental ^ fields, the following new phenomena take place. 
These fields are locally confined in the effective collective field theory, which emerges after 
integrating out the gauge fields. In addition to the old scalar density p x (4>), there is a 
flavor matrix density cr l J(<ft), which describes the \& condensate a x (X) = ^^(A — & x )^x 
with the color indexes summed, but the rest of indexes fixed. There are simply no objects 
in this collective field theory with \I/ and ^ at different points, in virtue of the gauge 
invariance. In terms of Wilson loops this local confinement corresponds to the zero area 
law W(C) = S 0tA (c)- 

Clearly, these $ fields cannot be quarks, they play different role. Namely, they break 
the local invariance of the scalar model. In the limit of Nf N, where the mixed 
model is solvable, this breaking is proportional to jf, unless there is a dynamical en- 
chancement. As we noted in [11|, such enhancement could, indeed, take place. The 
breaking combination of these densities satisfies the linear homogeneous integral equation, 
with the kernel depending on the old density. Under certain conditions, which generically 
involve the extra parameter adjustment, this linear integral equation would have finite 
solution. This is just the symmetry breaking we are looking for. In the limit of Nf N, 

3 At D — * 1 + one could test the solution, by perturbing the well-known free fermion Hamiltonian 
by the spatial interaction with the gauge field. This would be an interesting exercise, which might shed 
some light on the mysterious mechanism of the induction of the gauge theory. 
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when the model is solved classically, this breaking is spontaneous. 

There is the well known " no-go" theorem due to Elitzur, which forbids the spontaneous 
breaking of any finite local symmetry. The point is, for the finite symmetry group one 
could trust the local group integrations, which eliminate the nontrivial representations 
from the set of local observables. However, we are dealing with exceptional case, when 
these group integrals diverge in the large N limit. In fact, all our critical phenomena 
are hidden in the Itzykson-Zuber integral over unitary group. Much in the same way as 
in the Eguchi-Kawai model, the spatial divergencies are imitated by the divergencies of 
these group integrals. This is the key to understanding the large N induced models. 

When the physical, light quarks are added to the mixed model, they already can 
propagate in the vacuum. The classical Wilson loop would now obey the perimeter law, 
since the global Zn symmetry is also broken. Roughly speaking, the light quark screens 
its Zn charge by picking up \I/ from the vacuum, plus some number of pairs. The 
decrement of the perimeter law could be very large, though, if the symmetry breaking 
is small. So, the Wilson loops in the local limit would satisfy the area law, if the string 
tension would stay less than the mass cutoff, so that the area term would overcome the 
perimeter term. 

It is still unclear, how could one get the area law in the mixed model. Clearly, it is 
impossible at the classical level of effective field theory of two densities, so we need the 
quantum effects to be enchanced at large distances. 

In this paper we are studying the mixed model in detail. 



In the next Section 2. we reproduce the basic equations of the mixed model ||16|| , with 
some generalizations made and some bugs removed. 

In the next Section 3. we go to the local limit, in the same way, as in the scalar model, 
and find exact solution for the fermion density. In the scalar density, we find a Bose 
condensation of the eigenvalues at the origin, leading to the nonperturbative phenomena 
in the infrared region. 

In the Section 4. we study the wave equation for the scalar branch of the spectrum. 
The Bose condensation leads to new mass unit, which scales as a fractional power of 
-J- — > 0. We argue, that the mass spectrum is infinitely raising, as fractional power of the 
principal quantum number. 

In the Section 5. we discuss the remaining problems, and suggest the flavor instanton 
mechanism of the law. 

In the Appendix we compute the lowest scalar masses and the WKB asymptotics of 
the higher ones in the mixed model. 
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2 Basic equations of the mixed model 



The functional integral of the mixed model reads 



J] / D^exp-JVtr [U{$ x )\ f D^exp tr [M($ x )* x * x 
x j j 

JJ f DVt xy exp AT tr 



(2.1) 



<xy> 



$> x n xy <$> y n xy + n xy ^ y H xy ^ x + nl y y x H yx v y 



where ty x (^ x ) is N ®N'(N' (&N) matrix. The second dimension N' = N s Nf where Nf is 
the number of flavors, and N s = 2 Li^J is the number of spin components. The matrices 



H 



| (1 + 7 M ) act on the spin components. 



As compared to the previous paper [ lql , we introduced the direct interaction ^M($)$ 
between the scalar and spinor matter, instead of just the mass term \I/M\I/. This does not 
spoil the solvability of the model. We also removed an unnecessary factor of A" from this 
term. 



The one-link integral 



dVt exp A" tr \® x m 2 & + m 2 H l2 $ x + a^i# 2 i^ 



(2.2) 



depend, in virtue of the gauge invariance, only upon the densities pi, p%, &i, cr 2 , where 



Px(A) = — tr 6{\-$ x ), 



(ia\a x (X)\j(3) = - tr S(X - ^TW- 



(2.3) 



(2.4) 



with the spin indexes a, (3 and flavor indexes i, j fixed. We shall treat a as N' <8> N' matrix 
and use the notation tr' for the corresponding trace. 

Changing the variables from $ x , ty x to these local densities, we obtain the extra term 
in effective Lagrangean 

8S eff (a, p) = N J d\p(X) (-N' In p(A) + tr' In <x( A)) . (2.5) 
As a result, we get the following classical equations for these two densities 

(The factor p(X) in denominator of r\ was unadvertedly skipped in the previous paper |16j .) 

Next problem is how to compute the variations of the one-link integral / in these 
equations. In the same way, as in the scalar model, the Schwinger-Dyson identities are 



5 



used to derive the integral equations for these two variations.^ The first Schwinger-Dyson 
equation for the scalar density is the same as in the scalar model 



/i — A A — R(fj,) J 

but with the extra terms in the effective potential 

/N U'iv) N f N s ( d(v) a'(v)\ D-l mrrt/ . . . 

where 

V'(z) = I dfi-^-. (2.10) 

J Z — [A 

The second Schwinger-Dyson equation reads 

I Aii pWHy) , t, 1 (A - J?(^)) f dnjifj) 

J A — v J (a — R{is)) 2 + n z p 2 (v) J n z v — /i 

where 

G ° M = rW : ^ W = exp [I^T) arctan 1=550 J ■ (2 ' 12) 

The Schwinger-Dyson equations also yield the auxiliary condition 

dvr]{y). (2.13) 
These equations are exact in the large N limit, for arbitrary potentials U(X), M(X). 



3 Bose condensation 

For any particular lattice model, there would be a finite support of the eigenvalues, with 
the endpoint at the lattice cutoff scale. In the local limit, we are going to ignore these 
lattice artifacts, and extend the support to the whole real axis. The ultraviolet divergen- 
cies would be dimensionally regularized. As a result, there would be simple dispersion 
relations between the real and imaginary parts of the functions, such as V'(z), which are 
analytic in the upper half plane, with the symmetry 

V'(-z) = -V\z). (3.1) 

These dispersion relations read 

&V'(u + i0) = p 7r 2 pO = P / - -■ (3.2) 

J V — jl J V — jJ, 

4 We assumed that the classical solutions for 77(A), <r(A) are proportional to the unit N' ® N' matrix 
in virtue of symmetry. 
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Also, note, that in virtue of the symmetry of the scalar model, the density cr(A) is an 
even function of A as well, as the old density p(A). Therefore, rj(X) is even as well, so that 
the principle value integrals in ( [2.11| ) are the odd functions. The product of remaining 
factors in the v integral on the right of ( |2.11| ) are neither odd nor even; hence, only the 
odd part of this product can be left. 

With this property in mind, we can expand this product in powers of p and r(A) = 
-R(A) — A in the local limit, in the same way, as in the previous works. The first odd term 
would be a linear one, and we find 

P J du ^r = P J du J (,-0 )( A-0) P J ^—p> (3 - 3) 

Integrating over <fi and v on the right side, and using dispersion relations, we find a 
very simple result 

dpr)(p) 



•s 



A — p 

which leads to the puzzling equation 



(p(li) - p(X)) . (3.4) 



= p(A)p/^. (3.5) 
J A — p 

How could it be satisfied for finite rp. Recall that all our integrals are understood in 
the sense of analytic continuation of the powerlike integrals. For the powerlike Anzatz 
r](p) = p 1 we find 

= p(A)A 7 tan^ (3.6) 

This equation would hold for any even 7. Positive values of 7 correspond to nonsingular 
density rj which is dominated by the lattice scales. The singular solutions, with 7 = —2n 
are more interesting. With proper normalization this corresponds to 

dp^- oc iz~ 2n - p [ dp^- oc 3 (A + i0y 2n oc ^-^(A). (3.7) 
p — z J p — A 

Coming back to ( |3~5|) , we see that this solution holds for the density p(A) which 
vanishes faster than A 2n_1 at the origin. Note, that the auxiliary condition is also satisfied 
for any positive integer n. 

Now, let us substitute this solution into the first classical equation and find the first 
-J- correction to the scalar density. Simple algebra yields, up to irrelevant regular terms, 
at v — * 

N f N a ( p\v) a'{y)\ nN f N s l 

2ND\ p{y) a{v) ) ND v 1 ] 

This term can be compensated in R[y) by the last term, if the following 5-term would be 
added to the scalar density 

= w% s(x) - (3 - 9) 
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This term is positive, as it should be. This is the Bose-condensation of the eigenvalues. 

Note, however, that the sign changes at D < 1. In this case, the gap must develop, to 
avoid negative density. Or, to be more explicit, the pole term in R(u) at D < 1 can be 
compensated by the extra term in the density 

5p(\) = -p(X)<d(b 2 - A 2 ) (3.10) 

which effectively cuts the dispersion integral at ±6. The corresponding variation 

D — 1 r 5p(u) 2(D - l)b a+1 1 

p / dv^P^- -> —: — -r (3.11) 



D J X-u D(a + 1) A 

would cancel the pole term provided 

x = (a + l)nN f K 
2(1 - D)N 



(3.12) 



This is an interesting phenomenon. The transition at D = 1 corresponds to the 
interplay between the repulsion of the eigenvalues, caused by the usual Vandermonde 
determinant, and the attraction, caused by the gauge field via the Itzykson-Zuber deter- 
minant. Namely, the increase in the density at the origin raises the Vandermonde term 
in effective energy and lowers the Itzykson-Zuber term. At D > 1 the latter wins. As for 
the interaction with fermions, it effectively raises the energy, due to the Fermi repulsion. 
The net classical forces are balanced, therefore, by compensating this repulsion by the 
Itzykson-Zuber attraction. 

Let us now look at the local limit of the first classical equation |2] 

2r(A)p(A) = pj (3-13) 

Perturbing p in this equation (r(A) = R{\) — A is not perturbed), we find in the first order 

«W = P /*^. (3.14) 
J p — A 

This equation is satisfied for the Sp(X) oc <5(A) because r(0) = p(0) = p'(0) = 0. The 
vanishing first derivative p'(0) = is also needed, otherwise the right side would not 
vanish at A = 0. Vanishing of higher derivatives of p(X) at the origin is not needed here. 

The ^-function at the origin is probably just the result of the expansion. We expect 
this delta function to be smeared for finite jf. Actually, the shape of the density in the 
region A a+1 ~ -J- is beyond the perturbation theory. In this region, the higher orders of 
the -jf expansion are important, and we expect some scaling function 

V'(z) = z a f (^h) , m Jo ■ I,r • .... (3.15) 
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In other words, there is an infrared catastrophe in this model, taking place at small 
scalar fields, which is equivalent to small mass scales. Below the "confinement scale" 



/ Nf \ «Ti 

z ~ (^) (3.16) 

the quantum effects are significant for arbitrarily small ^f. Hopefully, the dominant terms 
are easier to sum up, than the QCD planar graphs. 



4 Scalar wave equation 

Let us study the wave equation for the fluctuations of the scalar density. The wave 
equation was derived in || and studied in JTT]. It reads (at z = A + zO) 



- ~ (P 2 + M 2 eff ) = D UV £F - %V + ^ D D l) \ v ( m ^ F ) , (4-1) 
where the analytic function 

= {-izT , (4-2) 



cos^f 



is the same as before. 



P is the Euclidean 4-momentum of the vacuum excitation, corresponding to plane 
wave fluctuations of p, 

Ap (A,x) = Tf e lPx p / dv \ , ' -. 4.3 

N J Ti 6 p{v) v — A 

The effective mass M 2 ^ gets the double pole term 

~ . . DU"(z) — 2D 2 — 2D + 2 a 2 a 2 
MZ f f{z) = D _ 1 + ^ = ^ + ^o + .-., (4.4) 

with 



2 2DnN f N s 

a = —, —r, (4-5) 

N(D - 1) ' K J 

which is the net result of the Bose condensation. The higher terms of the expansion are 
irrelevant, since the dominant z ~ -%= are small at Nt <C N. 

The boundary condition at A — > is that p(A)SjF(A + lO) = 0, to kill the 5(A) term, 
which comes from the factor A in the ratio of the real to imaginary parts of V. This 
boundary condition allows for the strong singularity ^J-(z) oc z w , uj > —a. However, with 
the 5-like perturbation in p, this is replaced by a weaker condition u > 0. 

Unfortunately, we cannot impose the boundary condition at z — 0, because of the 
abovementioned infrared catastrophe. All we can do at this point is to impose some 
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matching condition at z ~ zq, where the higher corrections to the density become im- 
portant. We could demand, that QJ 7 vanishes at zq, or take any other finite value. As 
we shall see shortly, the generic asymptotics for T is exponentially large, therefore, the 
actual matching value of T does not matter. Only the leading exponential term should 
cancel. 

Let us turn to the wave equation and consider the infinite sum of power terms as an 
Anzatz for T 

V sin^p} V ' 

where s = {0, 1} is the "A-parity" 

T{-z) = (-l)\F(z); Ap(-A, x) = (-l) s Ap(X, x). (4.7) 

Differentiating real and imaginary parts of V, J 7 , multiplying by and collecting the 
power terms, we find the following equation 

]T /(e) (X +a - l AC(e) + A e (P 2 + r ) + a 2 \^ 2 ) = 0, (4.8) 

e 

where 

C(e) = 2e cot {^y^j + 2 (^j^ + ^) tan ™ (4.9) 

The corresponding recurrent relation reads 

/(c) _ QP 2 + r )f(e + a - 1) + a 2 f(e + a + 1) ^ 

The leading power e satisfies the transcendental equation 

C(e o )=0, (4.11) 

which has infinitely many real solutions, approaching the integers at large eo- 
There are two scaling variables involved, 



(412) 

where the scaling function ^(u, i>) expands in double series 



ii _ 

g(u,v) = £ —r—— ( E»W0- (4-13) 

n =o sin (J(s + (a - l)n - e ) J ;= 

The expansion coefficients y?„(/) satisfy recurrent relation 

,n (1\ fn-ljl) +<fn-l(l- 1) 

^ (Z) = -C( eo -n(a-l)-2Z) (414) 
10 



with initial condition 

<p (l) = V (4.15) 

which is suitable for numerical or analytical computations. One can estimate, that tp n (l) 
decreases as ^, which is sufficient for the absolute convergence of the expansion. 

The spectrum is supposed to be defined from the boundary conditions at 

z ^ Zo =[a) • (416) 

This would yield the following mass spectrum 

-P' = Ml = r + A^pj\ L (s,e ). (4.17) 

with the critical index 

C = ^— ,0<C<1. (4.18) 
a + 1 

The scaling variables are finite, so that the expansion can be used to find the roots 
numerically. The results are presented in Appendix. They correspond to the lowest 
energy levels. 

The higher levels can be found analytically. In this limit the first variable tends to 
oo, and the second one tends to zero. Leaving the / = term , and using the WKB 
asymptotics for v? n (0), we find the following spectrum (see Appendix for the details) 

&(s,e ) ocL c ,L — oo (4.19) 



5 Discussion 

The mixed model looks much better than the original scalar one. The Bose condensation, 
caused by interaction with the fermions, leads to the nonperturbative phenomena in the 
infrared region. There is the new mass scale M, which goes to zero as a fractional power 
of the ratio of flavors to colors M 2 oc (Jjf^J ■ Beyond this scale one cannot neglect higher 
terms of the expansion. 

There are two sources of these terms. Some terms come from the nonlinearities of 
the classical equations for the density. Those are easy to account. The difficult terms 
come from the higher quantum corrections, starting from jfi. At 1 <C Nf <C N only 
the fluctuations of the fermion condensate are important, those of the scalar density are 
damped by the factors. It is a new interesting problem to sum up all the relevant 
infrared terms; at present we do not know how to do it. 

It is instructive to compare this situation with the one at large when we cannot 
solve the model, but know that it induces QCD. In this limit the induced gauge coupling 
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constant goes to zero as -j^-, which leads to exponentially small physical mass scale in the 
lattice units. The scalar field and its self interaction is irrelevant in this limit. 

When the ratio decreases, this scalar interaction becomes important. In absence 
of it, there would be a phase transition, at -d- = y when the induced self interaction of 
the hard gluons takes over the screening of the running gauge coupling by the matter. At 
smaller ^ this coupling would run to zero in the ultraviolet region. 

The induced QCD scenario corresponds to the situation, when strong enough scalar 
coupling prevents this phase transition, pushing it all the way to = 0. Then, at 
small -jf we would still be in the QCD phase. This possibility would display itself as 
the infrared catastrophe. At small enough mass scales the nonperturbative effects would 
become important. If our scaling solutions correspond to this phenomenon, then we could 
estimate, that effective gauge coupling goes to zero, as inverse logarithm of -jf 

Nf 

The singular behavior of the effective gauge coupling and the mass scale with ^ does 
not qiuite meet the old expectations. From the point of view of the usual RG analysis, 
this singularity could come about due to the new fixed point. The small perturbation 
from the fermions could drive the scalar coupling to this new fixed point, in which case 
the correlation length would, indeed, scale as some fractional power of ^f. The physical 
meaning of this mysterious phenomenon is yet to be understood. 

How are the gluons reproduced in terms of the effective field theory of the gauge 
invariant local composite fields er(A),p(A)? Clearly, the classical solutions for p,a do not 
describe the long range fluctuations of gauge field. This must be either the instantons of 
this effective theory, or the quantum fluctuations. At small -J- the quantum fluctuations 
are damped, so, perhaps, the flavor instantons could be responsible for gluons. 

Take the Wilson loop, for example. The space independent classical solution would 
simply yield the product of classical values of the link matrices, thus leading to perimeter 
law. As for the instanton, it potentially could yield the area law, as we know from the 
variety of models in lower dimensions. All it takes, is the distortion of the flavor matrix 
field cr x (\) in the region of the topology of the disk, bounded by the loop. Say, for the 
boundary conditions with the 27m rotation in one of the U(l) flavor subgroups along the 
loop, there must be instantons, with the classical action, proportional to the area, as it 
follows from the usual homological arguments. 

Leaving aside the Wilson loop, and the stringy mesons it is supposed to describe, let us 
turn to the Goldstone-like mesons, corresponding to collective excitations in our vacuum 
densities a, p. The p excitations were studied in this paper. The a excitations could be 
studied by the same method. Unlike the (simplest fixed points of the) scalar model || |TT[ . 
the mixed model has infinite mass spectrum, raising as fractional power of the principal 
quantum number. This is a good sign, supporting the conjecture of induction of QCD. 
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At this point it is too early to make definite conclusions. The euphoria of the first 
months of playing with the new toy has passed, but it would be equally silly to turn away 
from this model, just because it is not quite QCD. In fact, this is a whole new class of 
solvable models, which deserves serious study, as it teaches us nonperturbative Quantum 
Field Theory. As a by-product, we might induce QCD, or pass the D = 1 barrier in the 
string theory, or find something completely unexpected. 
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A Spectral equation 



We imposed the boundary condition SjF = at a particular point z = |zo, and com- 
puted the spectrum numerically for D = 4, a = 3.36901. The Taylor expansion for the 
logarithmic derivative ^paSlnjF in powers of x — (P 2 + To) with 50 digits of accuracy in 
coefficients was converted to the (9,9) and (10,10) Pade approximants, after which the 
roots of denominators were found. The accuracy Axj of the poles Xi of this function was 
estimated as the difference between the roots of denominators of these two approximants. 

Here are the values of e : {{x\, Axi}, . . .} for the first ten levels of e . Only real roots 
with < .01 are left. The empty list implies that no such roots were found in our 
approximation. 



2.05 
4.07 
6.09 
8.11 
10.1 
12.1 
14.1 
16.1 
18.1 
20.1 



{{-53.3, -0.0279}, {-36.6, -0.0000478}, {12., 1.59 lO" 11 }, {27.8, 1.07 10~ 7 }} 
{{-21.5, -6.14 10- 10 }, {31.8, 2.4 10~ 7 }} 
{{21.1, -1.29 10" 9 }} 
{} 

{{23.2, 5.86 10~ 9 }, {51.8, 0.00266}} 

{{10.7,-0.000486}} ^ ' ' 

{{2.83, -2.22 10" 10 }} 
{{-5.11, 2.66 10~ 8 }, {-3.71, 7.85 lO" 11 }, {8.66, -6.86 10~ 7 }, {17.9, -0.00165}} 
{{-2.54, -1.2 10' 11 }, {3.41, -5.18 HT 12 }, {20.5, -0.0000458}} 

{} 

The large roots cannot be computed by this method, as the convergence slows down. The 
WKB method can be applied in this limit. Let us set A — 1, and neglect r compared 
to P 2 . The asymptotic form of <p n (l) = <p (e — (a — l)n — 21) follows from the recurrent 
equation 

0(e) ~ r ( (a+ C ^( e) ) (~( a + 1 ) C "( £ ))^ T ; e = e -(a-l)n- 21. (A.2) 

The sum over I converges at I n, since the second argument of the scaling function 
is small at large | P 2 \ . The remaining sum over n is dominated by the pair of two complex 
conjugate saddle points in e plane, given by equation 

2 

(In 0(e))' = In (A.3) 

z 

At large complex e 

C"(e) - 2 tan (™) - 2z Sign (9fe) (A.4) 

which yield complex conjugate constants C, C' for the two saddle points e, e. As a result, 
we find 

a+l 

a 2 / P 2 \— ,, , 



15 



The logarithmic derivative can be estimated as follows 
d ft fee'*) 1 



<9P 2 \ / \ / (l _a)K(e^) a - 1 

where $ is the phase of the integrand at the saddle point, 



$ = ^fln0 + elnz — InP 2 ) (A.7) 

a — 1 



The spectrum is, therefore, quantized as follows, 



$ = tt(L + ^) (A.8) 



2 



After simple transformations we find 



27TI \ , „ . "+ 1 



(a + l)e I C I \ z a ~ x J V 2 

which yields the spectral equation in the text 
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